An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

可靠性设计优化的球空间设计法

可靠性设计优化（RBDO）在航空航天等领域已经被广泛应用，具有显著提高可靠性的效果；然而它包含两个迭代循环-总体优化循环和可靠性分析循环，双循环求解方法效率低下甚至难以承受。本文提出一种可以提高可靠性设计优化求解效率的方法-球空间设计法，该方法可以方便地表达搜索区域，最大可能失效点可以从极值点处精确获取；同时，该方法可以将总体优化以及可靠性分析同步进行。本文所提出的球空间法的有效性将通过2个算例进行测试，测试结果表明该方法不仅精度与同类方法持平，求解效率有明显提升，它将会对提高航空航天装备如卫星、火箭、飞行器等产品的可靠性提供更为有效的方法和技术支持。     

空间目标监视雷达信号频率的优选准则

随着在轨微小卫星及空间碎片数量不断增加，为保证在轨卫星的运行安全，需要对微小目标进行有效跟踪探测及稳定编目。雷达作为近地空间目标监视的主用设备，在对厘米级小目标进行探测跟踪时，其工作频率是决定其能否可靠探测关键因素，高频有利于探测小目标，并可获得相对稳定的RCS，但高频信号的波束小、搜索能力较低，实现难度大，也使得成本过高。为适于实际工程应用，需从效费比出发，对频率进行优选。本文从目标RCS的角度提出一种空间目标监视雷达信号频率的优选准则，给出该准则下周长波长比ka的取值范围。为了更加清楚地说明该优选准则的应用，给出了直径5 cm和10 cm的目标在该准则下信号频率的取值范围，并与国外工作做了比较，证明该准则可行。     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

一种柔性空间双臂机器人的协同控制和避障方法

针对柔性空间机械臂在轨服务应用需求，提出一种基于刚体运动与柔性振动相耦合的空间双臂机器人协同控制方法.首先引入空间位姿变量的概念，构造出面向协同控制目标的Jacobian矩阵，建立柔性空间机器人系统的刚柔耦合动力学模型，基于指定的最小距离得到其运动学逆解，并根据系统动量矩守恒关系及系统的Jacobian矩阵，并根据机械臂末端的运动速度,然后采用阻尼最小二乘法得出关节角度，使柔性空间机器人能够有效完成协同控制和空间避障任务，并基于RecurDyn V7R5软件环境验证算法的正确性.最后，基于SolidWorks和ADAMS虚拟样机建立柔性空间机器人系统的立体CAD模型，并结合空间在轨搬运任务进行模拟仿真，柔性空间机器人关节操作和运动轨迹的仿真结果图验证了本文算法的有效性.     

Injection of 40-kHz-modulated electron beam from the satellite: II. Excitation of electrostatic and whistler waves

Beam-plasma interaction effects are studied during the active space experiment with electron and Xe-ion beam injections in an ionospheric plasma. Permanent 40-kHz-modulated electron beam injection occurs simultaneously with a xenon-ion beam injected by the Hall-type plasma thruster operating in a square-pulse mode (100/50 s for a job/pause duration). The unusual behavior of the background charged particle fluxes and wave activity stimulated during the beam-plasma interaction have been registered by the scientific instruments onboard Intercosmos-25 station (IK-25) and Magion-3 subsatellite. The longitudinal and electromagnetic wave instabilities and their mutual relationship are considered in order to explain the observed effects. The excitation of electrostatic waves by the electron injection has been considered for different resonance conditions near the linear stability boundary. Beam-driven electromagnetic instability is responsible for the backward-propagating whistler waves excited via cyclotron resonance. Competition of these two beam instabilities is one of the subjects of the present study.     

期望状态序列导向的深空探测器规划修复方法

Aiming at the challenges caused by the persistence, concurrency and energy consumption of probe actions, a plan repair method of deep space probe based on the expected state sequence is proposed. In this method, the expected state sequence is formed of the expected effect of the unfinished action and the expected precondition of the unexecuted action in the pre designed plan, according to the execution status of the action. The expected state sequence is an ordered set of states with mixed logic and energy, providing subgoals for plan repair and also transforming the plan repair problem into the state transition path searching problem. During the search, the plan repair strategy with energy supply priority is proposed, which separates the logic repair from energy repair to reduce the difficulty of solving the problem. And this method enables the probe to recover from plan failure autonomously. Finally, the effectiveness and rationality of the proposed method are verified through simulation by taking the Mars Orbiter as an example.     

利用地月间空间站的载人登月飞行模式分析

为提高空间站利用率，降低载人登月任务成本，有效开发地月空间，研究了基于地月空间不同轨道空间站的载人登月飞行模式。首先对比直接往返登月飞行模式，对基于空间站的载人登月飞行模式进行任务分析，通过空间站将载人登月任务解耦为载人天地往返任务和登月任务两部分；其次通过轨道设计和稳定性分析提出考虑登月任务需求的地月间空间站可运行轨道和停泊点；最后建立一套飞行模式评价模型，从速度增量需求、飞行时间、空间环境、登月任务窗口、测控条件、交会对接技术难度、后续任务支持性和任务可靠性方面对6种不同位置空间站的登月飞行模式进行分析和定量评价。评价结果表明基于L2点Halo轨道空间站的载人登月飞行模式为更优飞行模式。     

密封舱内漂浮小球运动规律的数值模拟研究

"天宫二号"空间实验室开展了首次人机协同在轨维修技术试验,机械臂操作终端抓取漂浮小球试验为在轨试验任务之一。为研究漂浮小球与机械臂操作终端运动特性之间的匹配性,文章建立了在轨微重力环境下小球运动的物理模型与数学模型,采用动网格方法对不同小球的运动特性进行数值模拟,得到小球的运动规律,由此确定小球的设计参数范围,作为机械臂操作终端系统方案设计的依据。最终,数值模拟结果与"天宫二号"空间实验室在轨试验数据的对比初步验证了该计算结果的正确性。     

27-day variation in solar-terrestrial parameters: Global characteristics and an origin based approach of the signals

The Earth and the near interplanetary medium are affected by the Sun in different ways. Those processes generated in the Sun that induce perturbations into the Magnetosphere-Ionosphere system are called geoeffective processes and show a wide range of temporal variations, like the 11-year solar cycle (long term variations), the variation of ～27?days (recurrent variations), solar storms enduring for some days, particle acceleration events lasting for some hours, etc.In this article, the periodicity of ～27?days associated with the solar synodic rotation period is investigated. The work is mainly focused on studying the resulting 27-day periodic signal in the magnetic activity, by the analysis of the horizontal component of the magnetic field registered on a set of 103 magnetic observatories distributed around the world. For this a new method to isolate the periodicity of interest has been developed consisting of two main steps: the first one consists of removing the linear trend corresponding to every calendar year from the data series, and the second one of removing from the resulting series a smoothed version of it obtained by applying a 30-day moving average. The result at the end of this process is a data series in which all the signal with periods larger than 30?days are canceled.The most important characteristics observed in the resulting signals are two main amplitude modulations: the first and most prominent related to the 11-year solar cycle and the second one with a semiannual pattern. In addition, the amplitude of the signal shows a dependence on the geomagnetic latitude of the observatory with a significant discontinuity at approx. ±60°.The processing scheme was also applied to other parameters that are widely used to characterize the energy transfer from the Sun to the Earth: F10.7 and Mg II indices and the ionospheric vertical total electron content (vTEC) were considered for radiative interactions; and the solar wind velocity for the non-radiative interactions between the solar wind and the magnetosphere. The 27-day signal obtained in the magnetic activity was compared with the signals found in the other parameters resulting in a series of cross-correlations curves with maximum correlation between 3 and 5?days of delays for the radiative and between 0 and 1?days of delay for the non-radiative parameters. This result supports the idea that the physical process responsible for the 27-day signal in the magnetic activity is related to the solar wind and not to the solar electromagnetic radiation.